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Abstract 

Weak lensing leads to the non-Gaussian magnification distribution of 
standard candles at given redsliift z, p{n\z). In this paper, we give accurate and 
simple empirical fitting formulae of the weak lensing numerical simulation results 
with the generalized Dyer-Roeder prescription. The smoothness parameter a 
essentially represents the amount of matter that can cause magnification of 
a given source. Since matter distribution in our universe is inhomogeneous, 
we can think of our universe as a mosaic of cones centered on the observer, 
each with a different value of a. We define the direction dependent smoothness 
parameter a via the Dyer-Roeder equation; there is a unique mapping between 
a and the magnification of a source. We find that the distribution of a at given 
z, p{a.\z), is well described by a modified Gaussian distribution. For the same 
matter distribution, i.e., the same p{a\z), different values of f2m and can 
lead to very different magnification distributions. 

Our formulae can be conveniently used to calculate the weak lensing effects 
for observed Type la supernovae at arbitrary redshifts. 
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1. Introduction 

The use of standard candles is fundamental in observational cosmology. The distance- 
redshift relations for standard candles enable us to determine the basic cosmological 
parameters Hq (the Hubble constant), Q^n (the matter density of the Universe in units of the 
critical density pc = 3Hq/{87:G)), and (the density contribution from the cosmological 
constant in units of pc)- 

At present, the best candidates for standard candles are Type la supernovae (SNe 
la), because they can be calibrated to have very small intrinsic dispersions at cosmological 
distances ( Phillips 1993| , [Riess, Press, fc Kirshner 1995|) . Two independent groups of 



observers ( [Perlmutter et al. 1998| , [Riess et al. 1998D have demonstrated that SNe la 



can be potentially powerful tools for cosmology. Their current results ([Perlmutter et| 



al. 1999| , [Schmidt et al. 1998|) seem to indicate a low matter density universe, possibly 



with a sizable cosmological constant. It is clear that the observation of SNe la can 
potentially become a reliable probe of cosmology. However, there are important systematic 
uncertainties of SNe la as standard candles, in particular, evolution and gravitational 
lensing. The two groups have either assumed a smooth universe in their data analysis, or 
included lensing effects in rudimentary ways. Since we live in a clumpy universe, the effect 
of gravitational lensing must be taken into account adequately for the proper interpretation 
of SN data. 

A number of authors have considered various aspects of the gravitational lensing of 
SNe la ([Frieman 1997| , [Wambsganss et al. 1997] , [Kantowski 1998[ [Holz 1998| , [Holz fc Wald 



1998| , [Metcalf 1998| , [Porciani fc IVladau 1998[ ). A realistic calculation of the weak lensing 



effect of standard candles was conducted by Wambsganss et al. (1997), who computed 
the magnification distributions of standard candles using a N-body simulation which has 
a resolution on small scales that is of the order the size of a halo. However, since the 
magnification distributions depend on the cosmological model and redshift, the numerical 
results can not be directly used to compute the effect of weak lensing for an observed SN la 
at arbitrary redshift. 

In this paper we derive accurate empirical fitting formulae of the weak lensing numerical 
simulation results for the distribution of the magnification of standard candles due to weak 
lensing; these simple formulae can be used to account for the effect of weak lensing in 
the analysis of SN la data. In §2, we give analytical formulae for the angular diameter 
distance of a standard candle in terms of the smoothness parameter a and constants which 
depend on redshift and the cosmological parameters Qm and Q\. In §3, we generalize the 
interpretation of the angular diameter distance obtained in §2 by allowing the smoothness 
parameter a to be a direction dependent variable, the direction dependent smoothness 
parameter; we extract the distributions of a from the magnification distributions found by 
numerical simulations. In §4, we give analytical formulae for computing the magnification 
distributions at arbitrary redshifts. §5 contains a summary and discussions. 
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2. Angular diameter distance as function of the smoothness parameter a 



In a Hubble diagram of standard candles, one must use distance-redshift relations to 
make theoretical interpretations. The distance-redshift relations depend on the distribution 
of matter in the universe. In this section, we express the angular diameter distance to a 
standard candle in terms of the smoothness parameter, a, which is the mass-fraction of the 
matter in the universe smoothly distributed (Pyer fc Roeder 1973| ). 



In a smooth Friedmann-Robertson- Walker (FRW) universe, d = 1 in all beams; the 
metric is given by ds^ = dt'^ — a'^{t)[dr'^ /{I — kr'^) + r'^{d9'^ + sin^ 6* where a{t) is the 

cosmic scale factor, and k is the global curvature parameter {Q^. = 1 — — Q/^ 
The comoving distance r is given by ( [Weinberg 1972[ ) 



cm 



r[z) 



smn 



\nA^'^ / dz' 



nrJl + Z'f + ^^, + ^.{1+ Z') 



'1/2 



where "sinn" is defined as sinh if VL^ > 0, and sin if < 0. If flk = 0, the sinn and fl^s 
disappear from Eq.([^), leaving only the integral. The angular diameter distance is given by 
dA^z) = r{z)/{l + z), and the luminosity distance is given by d^i^z) = (1 + 

However, our universe is clumpy rather than smooth. According to the focusing 
theorem in gravitational lens theory, if there is any shear or matter along a beam connecting 
a source to an observer, the angular diameter distance of the source from the observer is 
smaller than that which would occur if the source were seen through an empty, shear-free 
cone, provided the affine parameter distance (defined such that its element equals the 
proper distance element at the observer) is the same and the beam has not gone through 
a caustic. An increase of shear or matter density along the beam decreases the angular 
diameter distance and consequently increases the observable fiux for given z. (Schneider, 
Ehlers, & Falco 1992) 

For given redshift z, if a mass-fraction a of the matter in the universe is smoothly 
distributed, the largest possible distance for light bundles which have not passed through a 
caustic is given by the solution to the following equation: 



d 

dz 



dPA 

dz 



Da{z = 0) = 0, 



+ -dinm{l + zfD, 

dDA 



0, 



dz 



z=0 



Ho' 



(2) 



where g{z) = (1 + zf^l + + VLx[{1 + z)-^ - 1] ( |Kantowski T998|) . The = form 
of Eq. (0) has been known as the Dyer- Roeder equation ( Pyer fc Roeder 1973| , Schneider et 
al. 19921) . 



The angular diameter distance for given smoothness parameter a and redshift z, 
Da (a, -z), can be obtained via the proper integration of Eq.(2). However, for our purposes. 
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it is useful to fit the angular diameter distance at given redsliift, DA{a\z), to a polynomial 
in a, with the coefficients dependent on z and cosmological parameters. For the ranges of a 
and z of interest, the angular diameter distance given by Eq.(E]) can be approximated as 



D 



A[a\z) 



where do = -Da (a = 0), and 



do{z) + a{z)a + b{z)a + c{z)a 



(3) 



a{z) 

b{z) 
c(z) 



-do + 3do.5 - 3di + di„ 



3 

2 [do 
—do + di — a — b. 



3a 

2do.5 + di- — 



(4) 



with do.5 = Da^ch = 0.5), di = DA{<y = 1), and (ii.5 
difference between Eq.(|^) and the solution to Eq.(0 
ranges of interest (see Fig.3). 



= -DaI^ = 1-5). Fig.l shows that the 
is much smaller than 0.1 % for the 



The constants do, do.5, di, and depend on the redshift z and the cosmological 
parameters and they are easily computed by integrating Eq.(0). Table 1 lists do, 
do.5, di, and di.5 for various cosmological models at various redshifts. 

Note that Eq.(|]) is usually used to describe a universe with a global smoothness 
parameter a, with < a < 1, since a is the fraction of matter in the universe which is 
smoothly distributed. However, Eq.(|^) is well defined for all positive values of a. Since 
a essentially represents the amount of matter that causes weak lensing of a given source, 
and matter distribution in our universe is inhomogeneous, we can think of our universe 
as a mosaic of cones centered on the observer, each with a different value of a. This 
reinterpretation of a implies that we have a > 1 in regions of the universe in which there 
are above average amounts of matter which can cause magnification of a source (see §3). 



3. The distribution of the direction dependent smoothness parameter a 

In the previous section, we showed that we can separate the dependence of the 
angular diameter distance on the mass-fraction of matter smoothly distributed [a) from 
its dependence on the redshift and cosmological parameters (see Eq.(3)). Unlike angular 
separations and flux densities, distances are not directly measurable. Let us generalize the 
angular diameter distance obtained in §2 by allowing the smoothness parameter a to be 
direction dependent, i.e., a property of the beam connecting the observer and the standard 
candle. In order to derive a unique mapping between the distribution in distances and the 
distribution in the direction dependent smoothness parameter for given redshift z, we define 
the direction dependent smoothness parameter a to be the solution of Eq.(|D (or Eq.(3)) 
for given distance Da{z). 
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Table 1: The constants do, do.^- 
models and redshifts. 



di, and (ii.5 from Eqs.@ and (|), for various cosmological 



(r2m,^A) 




z = 0.5 


z = 1.0 


z = 1.5 


z = 2.0 


z = 2.5 


z = 3.0 


z = 5.0 


(0.4, 0.6) 


do 


0.29088 


0.39373 


0.43849 


0.46109 


0.47378 


0.48150 


0.49400 




do.5 


0.28766 


0.37929 


0.40897 


0.41577 


0.41329 


0.40692 


0.37412 




di 


0.28445 


0.36517 


0.38069 


0.37329 


0.35782 


0.34000 


0.27469 




di.5 


0.28127 


0.35138 


0.35362 


0.33352 


0.30708 


0.28015 


0.19307 


(0.2, 0.8) 


do 


0.30893 


0.43237 


0.49007 


0.52031 


0.53764 


0.54831 


0.56578 




do.5 


0.30694 


0.42198 


0.46655 


0.48168 


0.48361 


0.47940 


0.44655 




di 


0.30496 


0.41175 


0.44374 


0.44485 


0.43304 


0.41610 


0.34456 




di.5 


0.30298 


0.40166 


0.42160 


0.40976 


0.38577 


0.35808 


0.25795 



(0.2, 0) 


do 


0.27246 


0.36342 


0.40377 


0.42485 


0.43710 


0.44480 


0.45790 




do.5 


0.27116 


0.35793 


0.39252 


0.40722 


0.41294 


0.41415 


0.40342 




di 


0.26987 


0.35248 


0.38147 


0.39006 


0.38963 


0.38488 


0.35335 




di.5 


0.26858 


0.34709 


0.37062 


0.37336 


0.36716 


0.35693 


0.30743 


(1,0) 


do 


0.25485 


0.32929 


0.35952 


0.37434 


0.38255 


0.38750 


0.39546 




do.5 


0.24973 


0.31077 


0.32577 


0.32602 


0.32096 


0.31401 


0.28527 




di 


0.24467 


0.29289 


0.29404 


0.28177 


0.26599 


0.25000 


0.19725 




di.5 


0.23968 


0.27565 


0.26424 


0.24134 


0.21710 


0.19455 


0.12797 



Note that in numerical simulations, as in the real universe, we can have two lines of 
sight with the same fraction of smoothly distributed matter, but different distances to a 
given redshift z, because weak lensing depends on where the matter is, as well as what 
fraction of the matter is smoothly distributed. Our definition of the direction dependent 
smoothness parameter a implies that it is no longer simply the fraction of matter smoothly 
distributed, it also contains information on where matter is distributed. We can interpretate 
a as the ratio of the effective density of matter smoothly distributed in the beam connecting 
the observer and the standard candle and the average matter density in the universe, with 
the effective matter density corresponding to a given amount of magnification. Two lines of 
sight with the same fraction of smoothly distributed matter but different distances would 
have different effective densities of smoothly distributed matter, thus different values of a. 

For given redshift z, we expect a distribution in the angular diameter distance Da{z) 
because the distribution of matter between redshift zero and redshift z is inhomogeneous. 
We have parametrized matter distribution with a, the direction dependent smoothness 
parameter; Eq.(3) then tells us how the angular diameter distance depends on the matter 
distribution for given redshift z. 

Since the direction dependent smoothness parameter a describes the distribution of 
matter in an arbitrary beam, it is a random variable for given redshift. The direction 
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dependent smoothness parameter a depends on the matter density in the beam connecting 
the observer and the source, as well as how the matter is distributed in the beam; for 
matter smoothly distributed throughout the beam, a < 1 in underdense beams, while a > 1 
in overdense beams. Note that here we do not consider the possibility that a significant 
fraction of matter can be in point masses in some beams. 

The matter density field is Gaussian on large scales and non-Gaussian on small scales, 
thus we parametrize the probability distribution of a in a form resembling the Gaussian 
distribution (see Eq.(5)). 

Wambsganss et al. have found numerically the distributions of the magnifications of 
standard candles at various redshifts, with z =0.5, 1, 1.5, 2, 2.5, 3, 5, for Vt^ = 0.4, 

= 0.6 ( IWambsganss et al. 1997] , [Wambsganss 1999| ). To extract the distribution of a, 
we use = [-0^(0; = 1)/-D^(q;)] (see Eq.(H)), where Da^ch) is given by Eq.(2) or Eq.(3). 
We find that the distribution of a at given redshift z is well described by 



p{a\z) = Cn 



exp 



a 



w a'^ J 



(5) 



where Cnorm, ctpeak, w, and q depend on z and are independent of a. Fig.2 shows Cnorm, 
otpeaki w, and q as functions of z; the points are extracted from the numerical p{fi\z), the 



solid curves are analytical fits given by 



^peaky^j 

w{z) 
q(z) 



10 



-2 



0.53239 + 2.79165 ( - j -2.42315 i^-j + 1.13844 i^- 



1.01350 - 1.07857 
0.06375 + 1.75355 



1 

5z 
1 

5z 



2.05019 



5z 



2.14520 



4.99383 I — + 5.95852 , 

.5zJ \5z 



5zJ 

1x3 



0.75045 + 1.85924 ( - 2.91830 (^0^ + 1.59266 ^ 



(6) 



is the normalization constant for given z. The parameter apeak{z) indicates 



Cnorm 

the average smoothness of the universe at redshift z, it increases with z and approaches 
C(peak{z) = 1 at 2: = 5; the parameter w{z) indicates the width of the distribution in the 
direction dependent smoothness parameter a, it decreases with z. The z dependences of 
C(peak{z) and w{z) are as expected because as we look back to earlier times, lines of sight 
become more filled in with matter, and the universe becomes smoother on the average. The 
parameter q{z) indicates the deviation of p{a\z) from Gaussianity (which corresponds to 
9 = 0). 

Fig.3 shows p{cx\z) for z = 0.5, 2, and 5. The solid line is derived from p{fi\z) 
found numerically by Wambsganss et al. (1997); the dotted line is given by Eq.(|]), with 
coefficients from Eq.(BD; the dot-dash line shows the difference between the solid curve and 
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the dotted curve; the dashed hne shows the Gaussian distribution related to Eq.(||), 



p'^{a\z) = Cn 



exp 



a 



(7) 



Note that it is difficult to see the dotted lines (our empirical ffiting formulae), because they 
are so close to the solid lines (the numerical results). This is not surprising since we have 
a 4-parameter fit to a smooth bell-like curve. Comparisons with Eq.(^ show how much 
p(a|z) deviates from the Gaussian distribution. 



For the same matter distribution, i.e., the same 'p{oi\z\ different values of fi„ 
can lead to very different magnification distributions (see §4). 



and f2/ 



4. The magnification distribution due to weak lensing 

At given redshift the magnification of a source can be expressed in terms of the 
apparent brightness of the source £(a), or in terms of the angular diameter distance to the 
source Da{oi): 



C{a) 
^~ C(a = l] 



DAia) 



(8) 



where C{a = 1) and -0^(0; = 1) are the flux of the source and angular diameter distance to 
the source in a completely smooth universe, and a is the direction dependent smoothness 
parameter (see Eq.(3) and §3). Since distances are not directly measurable, we should 
interpret Eq.(|D as defining a unique mapping between the magnification of a standard 
candle at redshift z and the direction dependent smoothness parameter a a.t z; a 
parametrizes the direction dependent matter distribution in a well-defined manner. 

The distribution in the magnification of standard candles placed at redshift z is 

. I X .da DA{a = l) ( dDA\~^ 

p{ld\z) = p{a\z) — = p{a\z) 



djj, 2/i^/2 \ da ) 

where the parameters a, 6, and c are given by Eq.(||); they depend on fi^, ^a, and z. 

Fig.4 shows p{^Ji\z) aX z = 0.5, 2, and 5 for the cosmological model fi^ = 0.4 and 
= 0.6. The solid line is the p(/i|^) found numerically by Wambsganss et al. (1997); the 
dotted line is given by Eq.(^, with p{a\z) given by Eq.(|^). Note the excellent agreement 
between our empirical fitting formulae and the numerical results. 

Fig.5 shows p{^Ji\z) for three cosmological models at three different redshifts: (a) 
z = 0.5, (b) z = 2, and (c) z = 5. The three cosmological models are: (1) flm = 1 and 
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Qa = 0; (2) = 0.2 and Q/^ = 0; (3) and Qm = 0.2 and Q/^ = 0.8. We have computed 

using Eq.@, with p{a\z) given by Eq.(^. Note that Fig.5 and Fig.4 have the same 
matter distribution (the same p{a\z)), but different cosmological parameters. 

Models with different cosmological parameters should lead to somewhat different 
matter distributions p{a\z). It would be interesting to compare numerical predictions for 
p{(y\z) from N-body simulations for different cosmological models. In the context of weak 
lensing of standard candles, we expect the cosmological parameter dependence to enter 
primarily through the magnification /i to direction dependent smoothness parameter a 
mapping at given z (the same a corresponds to very different n in different cosmologies). 



5. Summary and discussions 

We have derived accurate and simple empirical fitting formulae to the weak lensing 
numerical simulation results of Wambsganss et al. (1997). These empirical formulae can 
be conveniently used to compute the weak lensing effect of standard candles for various 
cosmological models. Our formulation is based on the unique mapping between the 
magnification of a source and the direction dependent smoothness parameter a; a is the 
ratio of the effective density of matter smoothly distributed in the beam connecting the 
observer and the source and the average density of the universe, with the effective matter 
density corresponding to a given amount of magnification. We find that the distribution of 
a is well described by a modified Gaussian distribution; this is interesting since the matter 
density field is Gaussian on large scale and non-Gaussian on small scales. 

We have derived empirical fitting formulae for p{a\z) (see §3) from the numerical 
magnification distributions, p{fi\z), found by Wambganss et al. (1997) for fl^ = 0.4, 
= 0.6. For the same matter distribution, i.e., the same p{a\z), different values of Qm and 
can lead to very different magnification distributions (see §4). It would be interesting to 
see how p{a\z) depends on the cosmological model ( |Wang 1999| ). 

Our empirical formulae can be used to calculate the weak lensing effects for observed 
Type la supernovae in general cosmologies and at arbitrary redshifts. At redshifts of a few, 
the dispersion in SN la luminosities due to weak lensing will become comparable or exceed 
the intrinsic dispersion of SNe la (|Wang 1998| ). The Next Generation Space Telescope 
(NGST) can detect SNe la at as high redshifts as they possibly exist; while there are 
theoretical uncertainties on the estimated SN la rate at high z, it is likely that the NGST 
will see quite a few SNe la at redshifts of a few ( ptockman et al. 1998|) . Complementing the 
NGST search, possible new ground based supernova pencil beam surveys can yield between 
dozens to hundreds of SNe la per 0.1 redshift interval up to at least z = 1.5 ( |Wang 1998| ). 
The systematic uncertainties of SNe la as standard candles will likely be well understood 
within the next decade. 
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We have used the numerical results by Wambsganss et al. (1997) in deriving our 
analytical formulae, which accurately describe the non-Gaussian magnification distributions 
of standard candles due to weak lensing. We note that Frieman (1997) gave analytical 
estimates of the magnification dispersions due to weak lensing which are of the same order 
of magnitude as the numerical results of Wambsganss et al. (1997), even though he did not 
consider the non-Gaussian nature of the magnification distribution. It would be useful to 
have numerical results for cosmologies other than the Qm = 0.4, ^2^=0. 6 model studied 
by Wambsganss et al. (1997). Holz & Wald (1998) studied a few cosmological models, 
including = 1, fl\=0; they computed their magnification distributions assuming that 
the mass in the universe is distributed in unclustered isothermal spheres. It would be 
interesting to compare how different assumptions affect the numerical results. 

Note that we have used the magnification distributions calculated by Wambsganss et 
al. (1997) using a N-body simulation which has a resolution on small scales that is of the 
order the size of a halo. These magnification distributions contain the information on how 
matter is distributed, including the clustering of galaxies. By generalizing the smoothness 
parameter a to a direction dependent variable, we have been able to describe the weak 
lensing of standard candles in a simple manner, leading to accurate empirical formulae 
which can be easily used to calculate the weak lensing effect of type la supernovae. The 
distributions of a also contain the information on how matter is distributed. The derivation 
of the distribution of a from the matter power spectrum should reveal how the measurement 
of p{<y\z) (via p{fi\z)) can probe the clustering of matter and structure formation in the 
universe ([Wang 1999|) . 
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Fig. 1. — The relative differences between the exact angular diameter distance and the 
analytical approximation given by Eq.(|^) as function of the smoothness parameter a, for 
z=0.5 (solid line), 1 (dotted line), 1.5 (short dashed line), 2 (long dashed line), 2.5 (dot - 
short dashed line), 3 (dot - long dashed line), and 5 (short dash - long dashed line). Note 
that a < 1 in underdense beams, while a > 1 in overdense beams. 

Fig. 2. — The coefficients Cnorm, o^peak, w, and q as functions of z; the points are numerical 
results, the solid curves are analytical fits given by Eq.(^). 

Fig. 3. — The distribution of the direction dependent smoothness parameter, p{a\z), for (a) 
z = 0.5, (b) z = 2, and (c) z = 5. The solid line is derived from p{fi\z) found numerically by 
Wambsganss et al. (1997); the dotted line is given by Eq.(^, with coefficients from Eq.(^); 
the dot-dash line shows the difference between the solid curve and the dotted curve; the 
dashed line shows the Gaussian distribution given by Eq.(|^). It is difficult to see the dotted 
lines because they are so close to the solid lines. Note that a < 1 in underdense beams, 
while a > 1 in overdense beams. 

Fig. 4. — The magnification distribution of standard candles, p{fi\z), for z = 0.5, 2, and 5. 
The solid line is the p{fj^\z) found numerically by Wambsganss et al. (1997); the dotted line 
is given by Eq.(^, with p{a\z) given by Eq.(H). It is difficult to see the dotted lines because 
they are so close to the solid lines. 

Fig. 5. — p{fJ^\z) at z = 0.5 (a), 2 (b), and 5 (c) for three cosmological models: (1) = 1 
and = 0; (2) n^n = 0.2 and = 0; (3) and = 0.2 and = 0.8. Note that Fig.5 
and Fig.4 have the same matter distribution (the same p{a\z)), but different cosmological 
parameters. 
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Fig. 1. — The relative differences between the exact angular diameter distance and the 
analytical approximation given by Eq.(|^) as function of the smoothness parameter a, for 
z=0.5 (solid line), 1 (dotted line), 1.5 (short dashed line), 2 (long dashed line), 2.5 (dot - 
short dashed line), 3 (dot - long dashed line), and 5 (short dash - long dashed line). Note 
that a < 1 in underdense beams, while a > 1 in overdense beams. 



-14- 




Fig. 2. — The coefficients Cnorm, ctpeak, and q as functions of z; the points are numerical 
results, the solid curves are analytical fits given by Eq.(^). (a) Cnorm- 




Fig. 2.— (b) apeak- 
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Fig. 2. — (c) w. 




Fig. 2.— (d) q. 
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Fig. 3. — The distribution of the direction dependent smoothness parameter, p{a\z). The 
sohd hne is derived from p{fi\z) found numerically by Wambsganss et al. (1997); the dotted 
line is given by Eq.(^), with coefficients from Eq.(^; the dot-dash line shows the difference 
between the solid curve and the dotted curve; the dashed line shows the Gaussian distribution 
given by Eq.(|^). It is difficult to see the dotted lines because they are so close to the solid 
lines. Note that a < 1 in underdense beams, while a > 1 in overdense beams, (a) z = 0.5. 




Fig. 3.— (b) z = 2. 




Fig. 3. — (c) z = 5. 




Fig. 4. — The magnification distribution of standard candles, p{^\z), for z = 0.5, 2, and 5. 
The sohd hne is the p{fi\z) found numerically by Wambsganss et al. (1997); the dotted line 
is given by Eq.@, with p{a\z) given by Eq.(|5|). It is difficult to see the dotted lines because 
they are so close to the solid lines. 
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Fig. 5. — p{lJ'\z) for three cosmofogical models: (1) = 1 and fl\ = 0; (2) fl^ = 0.2 and 
— 0; (3) and = 0.2 and Qa — 0.8. Note that Fig.5 and Fig.4 have the same matter 
distribution (the samep(Q;)), but different cosmofogical parameters, (a) z — 0.5. 
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Fig. 5.— (b) z = 2. 
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